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摘要 : 由 于 横断 山地 区 特殊 的 地 理 效应 ,第 四 纪 冰 期 气候 的 反复 变化 对 该 地 区 植物 的 地 理 分 布 和 居 群 遗传 
结构 产生 重要 影响 。 为 了 揭示 该 地 区 分 布 的 物种 的 谱系 地 理 结构 ， 检 测 了 该 地 区 特有 的 、 生 于 湿地 环境 的 
植物 海 仙 报 春 (Primula poissonii) 的 叶绿体 trnL-F 和 trnT-L 两 个 序列 变异 。 共 对 13 个 居 群 167 个 个 体 的 
两 个 叶绿体 片段 进行 测序 ， 共 发 现 11 个 单 倍 型 。 三 种 单 倍 型 被 共享 ， 其 他 单 倍 型 都 只 存在 于 单个 居 群 内 。 
分 子 变异 分 析 (AMOVA) 结果 表明 分 布 区 内 海 仙 报 春 的 遗传 变异 主要 存在 于 居 群 间 (99.08% ), 且 居 群 间 
的 遗传 分 化 很 高 (Gs,=0.916, Fop=0.99077, N=0.973), 有 着 极 显 著 的 谱系 地 理学 结构 (Nop >Gop, P< 
0.01) 和 较 低 的 居 群 间 平均 基因 流 (N, =0.08)。 结 合 般 套 分 支 分 析 (NCA) ， 本 文 的 研究 结果 推测 该 物种 
在 冰期 可 能 存在 多 个 避难 所 , 并 没有 表现 出 大 多 数 温带 植物 在 第 四 纪 冰 期 所 经 历 过 的 种 群集 体 扩张 和 迁移 
的 现象 ， 异域 片 段 化 被 认为 是 该 物种 现 有 单 倍 型 分 布 格局 形成 的 主要 原因 。 海 仙 报 春 的 谱系 地 理 结构 是 气 
候 的 反复 波动 和 横断 山特 殊 的 地 理 环 境 互相 作用 的 结果 。 
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Phylogeographical Study on Primula poissonii (Primulaceae ) 
from Hengduan Mountains * 


Min-Shu SONG'”, Ji-Pei YUE’, Hang SUN’, Zhi-Min LI'™ 
(1 School of Life Science, Yunnan Normal University, Kunming 650092, China; 2 Key Laboratory of Biodiversity 
and Biogeography, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650204, China) 


Abstract; The topological effect of the Hengduan Mountains (HM) during the Quaternary climatic oscillation must 
have important effects on spatial distribution and genetic structure of plants distributed there. In this study, we ex- 
amined the chloroplast trnL-F and trnT-L sequence variation of Primula poissonii , a relatively common alpine perennial 
endemic to this region. Eleven haplotypes were detected by analyzing 167 individuals from 13 populations of this spe- 
cies. Only three of them are commonly shared while the others are respectively fixed in the single population. An ana- 
lysis of molecular variance (AMOVA) for populations of P. poissonii showed that the genetic variation mainly resided 
among populations (99.08% ) , and the estimates of interpopulation differentiation were very high (G,,=0.916, F= 
0.99077, N,,=0.973). Further, a highly significant phylogeographic structure was present ( Ner>Gsr, P<0. 01) 
whereas the value of average gene flow was low (N,,=0.08) in the entire geographical range. Along with the nested 
clade analysis (NCA), our results indicate there were multiple refugia for this species during the glacial stages. We 
failed to detect the interglacial or postglacial range expansion in this species as commonly revealed for the other tem- 


perate plants. The nested clade analysis indicated that allopatric fragmentation was likely the major processes that 
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shaped the present spatial distribution of haplotypes in this species. The special phylogeographic structure of P. pois- 


sonii may have resulted from a combination of both climatic oscillation and complex topology of HM. 
y P. potogy 


Key words: Hengduan Mountains; Primula poissonii; Phylogeography; Chloroplast DNA; trnL-F; trnT-L 








横断 山 位 于 青藏 高 原 东 南部 ， 通 常 为 四 川 、 
云南 两 省 及 西藏 自治 区 东部 一 系列 南北 向 山脉 的 
总 称 ， 其 山岭 海拔 2 000 ~ 6 000 m， 岭 谷 高 差 一 
般 在 1000 ~2 000 m 以 上 ( 施 雅 风 等 , 1998)。 横 
断 山 地 区 地 形 复杂 ， 多 数 地 区 存在 纵向 平行 的 山 
脉 及 被 其 分 隔 而 成 的 河谷 ， 从 而 使 其 生境 强烈 斑 
块 化 。 该 地 区 的 是 全 球 生物 多 样 性 的 热点 地 区 之 
— (Boufford and van Dijk, 2004) , 植物 区 系 十 分 
丰富 ， 据 研究 与 统计 约 有 189 科 、1 500 属 和 8 
590 种 ， 其 中 2 783 种 为 该 地 区 特有 ( 李 锡 文 和 
李 捷 ，1993; Zhang 等 ，2009; Zhang and Sun, 
2010) 。 研 究 表明 ， 第 四 纪 时 期 这 一 地 区 经 历 了 
4 次 大 的 冰期 和 间 冰 其 作用， 每 次 冰期 和 间 冰 期 
的 作用 都 对 这 一 地 区 南北 成 分 融 汇 交 流产 生 了 深 
远 的 影响 ， 进 而 使 之 成 为 现代 温带 和 高 山 植物 区 
系 的 中 心 ( 刘 淑 珍 等 , 1986; Wu, 1988), iF 
多 重要 的 植物 类 群 尤其 是 高 山 植 物 如 杜 静 花 
(Rhododendron) 、 龙 胆 属 ( Gentiana) 虎 耳 草 
(Saxifraga) 、 风 毛 菊 属 ( Sausurea) 及 报 春 花 
(Primula) 等 世界 性 大 属 的 现代 分 布 中 心 (吴征 
Ri, 1979, 1980, 1987; 王 文 采 , 1992; 李 锡 文 和 
李 捷 , 1993 )。 尽 管 我 国 在 冰期 并 没有 像 欧 洲 、 
北极 或 者 北美 那样 发 生 大 规模 的 冰 关 ,但 是 在 青 
藏 高 原 地 区 这 种 气候 波动 影响 了 植物 的 分 布 范 
围 ， 使 其 随 冰 期 气候 的 波动 发 生 进 退 。 青 藏 高 原 
的 隆 升 和 第 四 纪 冰 川 造成 的 巨大 的 环境 变化 不 仅 
改变 了 植物 分 布 格局 ， 而 且 也 改变 了 物种 种 内 的 
遗传 结构 〈 王 献 溥 和 刘 玉 凯 ，1994; 孙 航 , 2002; 
Wang 等 , 2005; Liu 等 , 2006; Zhang 等 , 2005 ) 。 
同时 ， 吴 征 匀 (1980, 1987) 提出 横断 山地 区 是 
第 四 纪 冰 期 北 温带 植物 的 避难 所 ， 冰 期 后 成 为 现 
代 北 温带 植物 的 重要 起 源 地 和 辐射 地 之 一 。 
依托 分 子 生 物 学 技术 建立 起 来 的 谱系 地 理学 
(Phylogeography) 主要 是 通过 直接 检测 物种 不 同 
居 群 DNA 变异 和 谱系 分 析 来 探讨 该 物种 曾经 发 
生 的 种 群 扩 散 以 及 迁移 路 线 ， 从 而 进一步 推测 该 
物种 当前 分 布 格局 的 历史 成 因 ， 特 别 是 能 揭示 出 
由 于 受 第 四 纪 冰 期 与 间 冰 期 反复 交错 的 历史 事件 
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影响 。 利 用 这 种 手段 将 有 助 于 了 解 横 断 山 地 区 植 
物 分 布 格局 以 及 可 能 存在 的 冰期 避难 所 (Cheng 
等 ,2005; Ikeda and Setoguchi, 2007; Meng 等 ， 
2007; Wang 等 , 2009) DNA 测序 技术 是 目前 谱 
系 地 理学 研究 中 运用 最 广泛 的 技术 之 一 ， 而 植物 
谱系 地 理学 研究 主要 集中 研究 叶绿体 DNA 
(chloroplast DNA, cpDNA) 片段 ， 该 基因 组 具有 
TTE, FRR, Bie Ul, 、 结 构 简单 及 碱 基 
序列 重组 少 等 特点 ， 有 利于 遗传 分 析 。 叶 绿 体 基 
因 的 非 编码 序列 不 仅 适 用 于 分 子 系统 学 ， 同 样 适 
用 于 种 内 的 谱系 地 理学 研究 (Schaal 等 , 1998 ) ， 
并 且 由 于 cpDNA 的 单亲 遗传 使 得 其 基因 组 在 揭 
示 种 群 范围 变化 及 种 群 的 历史 动态 方面 非常 有 用 
(Newton 等 ，1999; Schaal 等 ,1998 ) ， 尤 其 有 利 
于 揭示 因为 环境 变化 或 冰期 影响 下 的 种 群 的 变化 
及 动态 (Bartish 等 , 2006; Petit 等 , 2003 ; Okaura 
and Harada, 2002; Fujii 等 ，2002; Newton 等 ， 
1999) 。 目 前 该 领域 的 研究 主要 集中 在 北美 及 欧 
YA, Beheregaray (2008) 的 综述 表明 ， 近 20 年 
的 谱系 地 理 研 究 文 献 有 77% 涉及 北美 洲 及 欧洲 
的 研究 ， 主 要 是 关于 第 四 纪 冰 期 避难 所 及 冰期 后 
种 群 扩张 的 研究 ， 而 中 国 西南 山区 及 喜马拉雅 地 
区 属于 研究 较为 薄弱 的 地 区 。 近 年 来 在 横断 山 及 
青藏 高 原 地 区 的 谱系 地 理学 研究 上 取得 了 不 错 的 
进展 ， 如 祁 连 圆 柏 (Juniperus przewalskii Kom, 
Cupressaceae) (Zhang 等, 2005) ,青海 云 杉 复合 
体 (Picea likiangensis Complex, Pinaceae) (Meng 
等 ,2007 ) 和 露营 乌 涉 (Aconitum gymnandrum , 
Ranunculaceae) (Wang “, 2009) 等 类 群 的 研究 
揭示 了 中 国 -喜马拉雅 地 区 的 末次 冰期 存在 多 个 
避难 所 ， 这 些 类 群 通过 改变 它们 的 分 布 区 域 来 啊 
应 第 四 纪 气 候 的 波动 。 而 对 条 纹 狭 蕊 龙 胆 (Meta- 
gentiana striata, Gentianaceae ) (Chen 等 ，2008 ) ， 
偏 花 报 春 (Primula sencudiflora, Primulaceae ) 
(Wang 等 , 2008) 等 草本 类 和 群 的 研究 表明 ， 横 断 
山地 区 是 第 四 纪 冰 期 的 避难 所 。 

本 研究 所 用 的 材料 海 仙 报 春 (Primula poissonii 
Franch. ) 隶属 于 报 春花 属 灯 台 报 春 组 (sect. Pro- 




































































1 期 宋 敏 舒 等 : 横断 山地 区 海 仙 报 春 的 谱系 地 理学 研 


ot 





93 





liferae) ， 多 年 生 草 本 ， 生 于 海拔 2 100 ~3 500 m 
的 山坡 草地 湿润 处 和 水 边 ， 该 种 植物 主要 在 横断 
山区 的 南部 和 东部 (四 川西 部 和 云南 西北 部 ) 分 
布 ， 是 典型 的 异 花 授粉 植物 。 由 于 草本 植物 较 短 
的 生活 周期 以 及 对 生活 环境 变化 的 敏感 性 ， 加 上 
该 种 特殊 的 湿地 生境 ， 所 以 十 分 适合 作为 横断 山 
区 的 一 个 模式 种 来 研究 高 山 植 物 的 谱系 地 理学 。 
在 本 研究 中 ， 我 们 选取 了 cpDNA 的 两 个 非 编 码 
片段 trnL-F 及 trnT-L 对 海 仙 报 春 13 个 居 群 167 
个 个 体 测 序 ， 进 行 谱系 地 理学 研究 ， 主 要 目的 是 
通过 初步 取样 的 研究 ， 揭 示 该 植物 单 倍 型 的 地 理 
分 布 式 样 、 居 群 遗传 结构 以 及 冰期 时 的 避难 所 及 
冰期 后 的 迁移 路 线 。 




















1 材料 和 方法 
1.1 实验 材料 

本 研究 材料 采 自 四 川西 部 及 云南 西北 部 ， 按 居 群 采 
样 标准 对 海 仙 报 春 所 处 地 理 分 布 范 围 进 行 取样 ， 即 每 个 
居 群 随机 采样 7 ~ 15 个 个 体 ， 每 个 居 群 内 采样 个 体 之 间 
至 少 间隔 10 m， 共 取得 海 仙 报 春 13 个 居 群 167 个 个 体 
(#1, 图 1)。 野 外 采集 生长 良好 ， 无 病 斑 的 幼 叶 ,用 
硅胶 快速 干燥 后 带 回 实 验 室 。 和 凭证 标本 存 于 中 国 科学 院 
昆明 植物 研究 所 标本 馆 (KUN) 。 选 取 同 组 的 亲缘 种 玉 
山 灯 台 报 春 (Primula miyabeana) 作为 外 类 群 。 










































































1.2 DNA 提取 和 PCR 扩 增 


采用 改进 的 CTAB 法 (Doyle and Doyle, 1987) 从 硅 
胶 干 燥 的 叶片 中 提取 总 DNA, 对 cpDNA 的 trnL-F 片段 和 
trnT-L 片段 进行 PCR 扩 增 。PCR 反应 体系 在 20 pL 体积 
中 进行 ， 包括 模板 DNA 10 ~40 ng, 2 kL 10xPCR 反应 组 
冲 液 ，1.6 uL MgCl, (25 mmol +L"), 1.6 pL dNTP 混合 
Wi (10 mmol - L`), IER SIH 1 uL (10 pmol +L"), 
Taq DNA 聚合 酶 0.5 单位 。 反 应 程序 如 下 : 56 94°C 条 件 
下 变性 4min， 而 后 94% 1min, 57C 1 min, 72°C 1 min, 
循环 34 次 ， 最 后 在 72% 条 件 下 延伸 7 min。 扩 增 得 到 的 
PCR 产物 在 1% 的 琼脂 糖 胶 上 电泳 ,在 ABI 3700 DNA 
analysis system 测序 仪 上 测序 。 
1.3 数据 分 析 

对 167 个 个 体 的 PCR 扩 增 片段 的 全 序列 测序 ， 用 
CLUSTAL_X 1.81 (Thompson 等 , 1997) 软件 对 测定 的 序 
列 进 行 对 位 排列 ， 并 加 以 手工 校对 。 用 DnaSP 4.0 (Ro- 
zas 等 , 2003) 程序 确定 cpDNA 单 倍 型 。 根 据 每 个 居 群 内 
不 同 单 倍 型 的 个 体 数 目 和 总 的 个 体 数 来 计算 单 倍 型 频 
率 ， 然 后 基于 每 个 居 群 的 单 倍 型 组 成 、 不 同 单 倍 型 在 每 
个 居 群 的 个 体 数目 和 单 倍 型 序列 变异 位 点 ， 分 别 计算 它 
们 的 遗传 多 样 性 大 和 H, (分 别 为 居 群 内 和 居 群 间 平 均 
遗传 多 样 性 ) 以 及 居 群 遗传 分 化 Gs (Raymond and 
Rousset, 1995) 和 Nw 值 ， 依 照 Pons 和 Petit (1996) 所 描 
述 的 方法 ， 用 PERMUT 软件 计算 。 居 群 间 的 Ns; 值 用 
1000 次 置换 进行 置换 检验 。 










































































R1 海 仙 报 春 13 个 居 群 材料 来 源 及 叶绿体 DNA 单 倍 型 (H1-H11) 的 遗传 多 样 性 (Hd) 组 成 和 频率 


Tablel Origin of materials, measuremnets of haplotypes diversity (Hd) and composition/frequency from 


cpDNA haplotypes (H1-H11) in 13 populations of P. poissonii 






























































编号 样本 数 纬度 海拔 单 倍 型 数目 单 倍 型 多 BERL 
编号 eE 采集 地 . . ne ear 
No. Population Locality Sample Longitude Latitude Altitude Haplotype 样 性 指数 样 性 指数 
size (N) (m) number Hd T 
1 JL 川 九龙 Jiulong, Sichuan 13 101°23’ 29°50’ 2695 H5 (13) 0 0 
2 MN pu )I] 2% Mianning, Sichuan 15 102°12’ 28°39’ 3357 H8 (15) 0 0 
3 ML-A ”四 川 森 里 Muli, Sichuan 7 101°10’ 28°04’ 2985 H5 (7) 0 0 
4 ML-B 四川 森 里 Muli, Sichuan 8 101°18’ 28°00’ 3100 H7 (8) 0 0 
5 WS 云南 全 上 Wengshang, Yunnan 15 99°45’ 28?15 3675 H9 (15) 0 0 
6 ZD 云南 中 旬 Zhongdian, Yunnan 15 99°52’ 28°07’ 3460 H1 (15) 0 0 
7 WX 云南 维 西 Weixi, Yunnan 15 99°16’ 27°19’ 2897 H10(14),H11(1) 0.133(0.11) 0.0001 
8 ”XZD-A 云南 小 中 徊 Xiaozhongdian, Yunnan 13 99°47’ 27°34’ 3256 H1 (13) 0 0 
9 ”XZD-B 云南 小 中 旬 Xiaozhongdian, Yunnan 15 99°46’ 27°35’ 3367 HI10(6),H1(9) 0.514(0.06) 0.00039 
10 EY 云南 洱 源 Eryuan, Yunnan 15 99°59’ 26°20’ 2301 H1 (15) 0 0 
11 LJ 云南 丽江 Lijiang, Yannan 14 100°11’ 26°59’ 2943 H6 (14) 0 0 
12 KM 云南 昆明 Kunming, Yunnan 8 102°37’ 25°03’ 1925 H2(7) ,H3(1) 0.250(0.18) 0.00038 
13 JD 云南 景 东 Jingdong，Yunnan 15 101°01’ 24°32’ 3349 H4 (15) 0 0 
Total 167 0.791(0.02) 0.00305 
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图 1 海 仙 报 春 居 群 采样 点 及 叶绿体 DNA 单 倍 型 分 布 图 ， 居 群 编号 (1 ~13) 与 表 1 的 编号 一 致 
Fig.1 The sampled locations and distribution of cpDNA haplotypes within and among populations of P. poissonii. 


All of the population codes from 1 to 13 are accordant with Table 1 





应 用 ARLEQUIN 软件 包 version 3.01 (Excoffier 等 , 2006) 后 从 采样 地 点 之 间 的 地 理 距离 计算 Clade distance (De) 
中 的 AMOVA (Analysis of Molecular Variance) 分 析 方 法 和 Nested clade distance (Dn) o De 是 同一 单 倍 型 (分 支 ) 
分 别 检测 材料 在 分 布 区 内 居 群 内 和 居 群 间 的 遗传 变异 水 的 个 体 与 此 单 倍 型 分 布地 理 中 心 的 平均 地 理 距 离 ; Dn 
平 ， 以 及 单 倍 型 分 布 的 Fs 评价 (Weir and Cockerham, 1984) , 是 同一 单 倍 型 (分 支 ) 的 个 体 与 上 一 层 分 支 地 理 中 心 的 
以 进一步 分 析 居 群 的 分 化 程度 (1000 次 置换 检验 ) 。 平均 地 理 距 离 。 在 内 套 分 支 范 围 内 没有 地 理 联 系 的 零 假 

应 用 PAUP” 4.0b10 (Swofford, 2003) 软件 用 最 大 简 设 用 Geodis 2.5 (Posada 等 ，2000) 软件 进行 分 析 。 通 过 
约 性 分 析 法 (maximum parsimony, MP) 进行 单 倍 型 系统 比较 末端 分 支 (tip clade， 位 于 网 状 进 化 树 边 缘 、 只 后 

之 育 分 析 ， 以 亲缘 种 玉山 灯台 报 春 (P. miyabeana) VE 一 个 单 倍 型 相连 的 分 支 ， 简 称 T) 与 内 部 分 支 (interior 
为 外 类 群 ， 空 位 (gap) 作为 缺失 处 理 ， 采 用 局 发 式 搜 clade， 位 于 网 状 进化 树 内 部 、 与 多 个 单 倍 型 相连 的 分 
索 ，1000 次 随机 加 入 序列 ，TBR 枝 长 交换 ， 每 步 保 留 10 SE, 简称 1) 的 De 和 Dn 观测 值 、 或 在 从 套 分 支 范 围 内 
棵 树 。 得 到 的 一 致 性 系统 树 分 支 的 可 靠 性 使 用 “ 靴 带 ” 用 1000 次 Monte Carlo 置换 的 零 分 配对 采样 地 点 的 内 部 - 
(Bootstrap) 分 析 ， 用 1000 次 重复 来 检验 单 倍 型 各 分 文 末端 分 支 (interior-tip clade, I-T) 的 De 和 Dn 比较 时 进 
的 支持 率 值 。 行 检 测 。 如 果 内 部 与 末端 的 De 和 Dn 差异 不 显著 ， 零 假 

种 内 单 倍 型 间 骨 套 分 支 分 析 (NCA) 参考 Templeton 设 不 能 被 推翻 ， 表 明 居 群 之 间 没 有 地 理 相关 性 ;， 如果 内 
等 (1992) 的 方法 ,利用 TCS version 1.2 (Clement 等 ， 部 与 末端 的 De 和 Dn 差异 显著 ， 零 假设 被 推翻 ， 表 明 居 
2000) 软件 进行 构建 。 这 种 方法 应 用 漳 祖 理论 ( Coales- 群 之 间 存 在 地 理 相 关 性 ， 则 需要 进一步 推断 居 群 之 间 的 
cent theory) 来 推测 单 倍 型 间 的 亲缘 关系 。 当 网 状 图 中 间 谱系 地 理 关系 。Templeton (2004) 的 检索 表 用 于 推断 单 
出 现 了 闭合 环 状 结构 时 ， 我 们 使 用 基于 溯 祖 理论 的 一 种 信 型 变异 之 间 的 可 能 的 时 空 联系 。 
方法 和 原则 进行 解释 (Crandall and Templeton, 1993; 物种 水 平 上 的 居 群 间 的 平均 基因 流 CN, ) 、 单 倍 型 多 样 
Templeton and Sing, 1993 ) 。 单 倍 型 间 的 地 理 相关 性 首先 性 指数 (Hd), AERE FENES Z (ar) (Nei, 1987) 以 及 
] Templeton 和 Sing (1993) 在 舱 套 范围 内 对 采样 地 点 的 Tajima’s D (Tajima, 1989), Fu 和 Li's D* (Fu and Li, 
EARAY SCA EFT REST BEE RM LY AE ETF, YAR 1993) 两 种 无 限 突变 位 点 模型 的 中 性 检验 方法 检测 及 失 配 分 
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析 ， 又 称 歧 点 分 布 (Mismatch distribution) 分 析 (Rogers 
and Harpending, 1992) 都 在 DnaSP 4.0 程序 中 完成 。 





2 结果 
2.1 序列 变异 

对 13 个 居 群 167 个 个 体 进行 了 叶绿体 非 编 
码 DNA trnL-F 和 rn T-L 片段 的 全 序列 测定 ， 通 
过 排序 发 现 ，tirnL-F 片段 长 度 为 940 ~ 946 bp, 
trnT-L 片段 长 度 为 448 ~ 460 bp， 两 片段 联合 长 
JE 1406 bp， 共 检测 到 11 个 单 倍 型 (H1-H11 ) 。 
这 11 种 单 倍 型 变异 位 点 见 表 2。 
2.2 ”多样 性 指数 、 遗 传 结构 分 析 及 失 配 分 析 

在 海 仙 报 春 的 13 个 居 群 总 的 单 倍 型 多 样 性 
SŽ Hd 为 0.791 (0. 023 ) ， 核 苷 酸 多 样 性 指数 
m 为 0.00305。 其 中 单 倍 型 Hl 的 频率 最 高 ， 有 
52 个 个 体 拥有 此 单 倍 型 。 每 个 居 群 的 遗传 多 样 
性 (Hd) 、 单 倍 型 数量 、 单 倍 型 组 成 及 频率 见 表 
1， 单 倍 型 的 地 理 分 布 情况 见 图 1 。 

通过 PERMUT 程序 计算 得 出 海 仙 报 春 居 群 
内 平均 遗传 多 样 性 H, 值 、 总 的 遗传 多 样 性 H, 
值 、 居 群 间 遗传 分 化 Gs 和 Nsy 值 分 别 为 :0. 078 
(0.0445) 0.925 (0.0431)、0.916 (0.0500) 和 
0.973 (0.0172)。 使 用 U- 统 计 方 法 对 海 仙 报 春 整 
个 分 布 区 单 倍 型 变异 的 地 理 结 构 进 行 检验 后 发 现 
Ns 极 显 著 的 大 于 Gy, (U=1.95，P<0.01) ， 表 明 
在 海 仙 报 春分 布 范围 内 ， 居 群 间 的 遗传 分 化 水 平 



































表 2 





很 高 (Gop =0.916)， 并 且 有 着 明显 的 谱系 地 理 
学 结构 出 现 。 分 子 变异 分 析 (AMOVA) 的 结果 
表明 ， 居 群 间 的 遗传 变异 为 99.08% ， 而 居 群 内 
的 遗传 变异 为 0.92% ，F,,=0.99077 (P<0.001) 
( 表 3) ， 进 一 步 揭示 了 海 仙 报 春 的 遗传 变异 主要 
存在 于 居 群 间 ， 而 且 具 有 较 高 的 居 群 分 化 水 平 。 
通过 DnaSP 4.0 程序 计算 得 出 居 群 间 的 平均 基因 
TA 0.08, 说明 居 群 间 基因 交流 不 频繁 。 在 种 
的 水 平 上 对 叶绿体 DNA 进行 的 失 配 分 析 ， 得 到 
的 失 配 曲线 为 多 峰 曲 线 (图 2) ， 背 离 了 对 快速 
扩张 模型 的 假设 ， 这 个 结果 在 两 种 基于 无 限 位 点 
0.25 9 
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Pairwise differences 


图 2 海 仙 报 春 所 有 个 体 叶 绿 体 DNA 两 个 片段 
联合 数据 的 失 配 分 析 

Fig.2 Mismatch distribution for sequences data of two cpDNA 

fragment from all individuals of P. poissonii. The thin represents 


the expectation; The dotted line represents the observed value 


海 仙 报 春 2 个 时 绿 体 片段 11 个 单 倍 型 间 的 序列 变异 位 点 


Table 2 Variable sites of the aligned sequences of two chloroplast DNA fragments in 11 haplotypes of P. poissonii 


i=} 


变异 位 点 Variable Sites 








单 信 型 1 1 1 2 2 2 2 234445777 8 8 9 1 1 1 1 1 1 
i 6 7 0 4 5 7 9 1 6 7 8 5 0 4 7 0 2 2 2 1 1 1 2 2 2 
6 4 5 2 2 1 5 5 0 5 8 7 2 4 2 5 7 8 7 1 6 8 1 4 4 
7 0 7 9 2 3 
H1 C G - C CTT- - CC C Ay - G - A AAA - G- A A A 
H2 T Cc - T C C T - - C G C Å, As A - A, A A Ay As A - Ap A A 
HB T C - € C CT- ~ C G A Ay As A - A, A A A, AS A - AS A A 
H4 T C - C C C T - - C G C Ay As A Ae A, A A A, å; A A, å; A A 
H5 c c - CTO T = = C A C Ay = G = = T AA = C= = TA 
H6 Cc G - C CTT -4 C G C A, - G - A A A A - G - AD A A 
H7 To ay © T & T = = © A € ee & e+ T Ce ve S&H + T €C 
H8 C Cc - C C CT- = T CGC C Ay å; A - A, A A Ay AS A - A, A A 
H9 C G - C C TT - - C GCA - G- A C A A, - G - & C A 
HI0 C G - C C T AA, - C G CA - G - A, A A A - G - Å A A 
H11 T G - C C TA- = CC CA = G - A, A A Å, - G - Å A A 
A,, TATTGG; A,, TTCTATGAA; A, , GAAGA; A, , AGTTTATIT; A;, TATTTA; As, AAACT; A, CTITTTGA; —, delete 
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表 3 所 有 单 倍 型 AMOVA 分 析 
Table 3 Analyses of molecular variance (AMOVA) for all haplotypes 











变异 来 源 Source of variation 自由 度 df 方差 SS 变异 成 分 VC 变异 比例 Variatioin (% ) 
居 群 间 Among population 12 1471.771 9. 58204 99.08 

居 群 内 Within population 154 13.750 0. 08929 0.92 

总 计 Total 166 1485.521 9. 67132 

国定 指数 Fixation index F op 0. 99077 ( P<0. 001) 





非 重组 模型 的 中 性 检验 Tajima’s D=0.94711 (P> 
0.10), Fu and Li's D=1. 10072 (P>0. 10) 中 得 
到 了 支持 ， 中 性 检验 结果 均 为 正 值 显示 我 们 所 采 
样 的 海 仙 报 春 种 群 由 于 长 期 的 地 理 隔离 增加 了 居 
群 间 的 分 歧 。 
2.3 各 单 倍 型 的 之 间 的 系统 发 育 和 艇 套 分 析 
(NCA) 

用 PAUP 软件 分 析 11 个 单 倍 型 ， 以 玉山 灯 
台 报 春 为 外 类 群 进行 系统 发 育 分 析 。 最 大 简约 性 
(MP) 分 析 ， FIERA 26 个 变异 位 点 ， 其 
中 10 个 信息 位 点 ， 得 到 一 棵 严格 一 致 性 树 (图 
3) 。 从 严格 一 致 性 树 看 出 ， 海 仙 报 春 13 个 居 群 
的 11 个 单 倍 型 分 化 为 3 个 主要 分 支 : 分 支 1 包 
括 单 倍 型 HI 、H6、H9 、HI10 及 HI1， 分 别 固定 
在 云南 西北 部 的 ZD, WS, WX, XZD, EY 及 
LJ; 分 别 固定 在 ML 两 个 居 群 和 JL 居 群 的 单 倍 
型 H5 和 H7 以 较 高 的 支持 率 (97%) RAM 
I, 分 支 亚 由 单 倍 型 H2、H3、H4、H8 组 成 ， 
分 别 固定 在 KM、JD 及 四 川西 部 的 MN 居 群 。 用 
TCS1.2 软件 构建 了 单 倍 型 间 的 最 小 生成 网 状 图 
并 进行 了 分 层次 的 般 套 进化 分 支 分 析 (图 4)， 
其 分 布 格局 与 最 简约 树 基本 一 致 。11 个 单 倍 型 
组 成 了 7 个 一 级 分 支 (1-step clades), 5 个 二 级 
分 支 (2-step clades), 3 个 三 级 分 支 (3-step 
clades), Clade3-1 包括 2 个 二 级 分 支 和 2 个 一 级 
分 支 ， 固定 在 ML 和 JL 两 地 的 3 个 居 群 中 ; 
Clade2-4 包括 2 个 一 级 分 支 Cladel -4 和 Cladel 
-5， 由 云南 西北 部 的 ZD、WX、XZD-A、XZD- 
B, EY, LJ 6 个 居 群 拥有 的 单 倍 型 组 成 ，Clade2 
-5 的 上 级 分 支 为 Clade3-3， 也 包括 2 个 一 级 分 
支 Cladel -6 和 Cladel -7， 分 别 由 KM、JD 及 
MN 所 拥有 的 单 倍 型 组 成 。 检 索 表 分 析 推 测 得 到 
的 各 个 支 系 单 倍 型 地 理 分 布 格局 的 一 系列 历史 成 
因 见 表 4: 整个 支 系 和 各 个 支 系 的 形成 原因 均 为 


















































异域 片段 化 (Allopatric fragmentation ) 。 


3 Wit 

从 图 3 可 以 看 出 11 种 单 倍 型 可 以 分 为 地 理 
范围 各 不 相同 的 3 SOP SC, Clade-I 分 布 于 云南 
西北 部 ， 其 中 单 倍 型 HI 被 相 邻 的 4 个 居 群 (中 
旬 、 小 中 旬 -A、 小 中 旬 -B 及 洱 源 ) 分 享 ， 网 状 
进化 树 支 持 其 可 能 为 这 个 分 支 的 祖先 单 倍 型 ， 这 
一 分 支 的 其 他 单 倍 型 都 通过 较 小 的 突变 从 这 个 祖 
先 单 倍 型 衍生 出 来 ， 并 且 固 定 在 不 同 的 地 区 ; 
Clade-I 仅仅 分 布 于 四 川西 部 的 木 里 -A、 木 里 -B 
及 九龙 居 群 ， 从 网 状 进 化 树 的 结果 看 来 ， 这 个 分 
支 经 过 了 较 多 的 突变 从 Clade- [中 分 化 出 来 ， 表 

















Primula miyabeana 


H2 昆明 








H3 昆明 
Clade-II 
H4 景 东 


H8 F 





HS 九龙 、 木 里 _A 





Clade-II H7 KELB 











H10 维 西 、 小 中 多 _-B 











H11 小 中 旬 _-A 

H1 中 多、 小 中 旬 _A、 
Clade-I 河源 
H6 丽江 
Ho E 


图 3 海 仙 报 春 11 种 单 倍 型 (HI ~ HI1) 的 1 棵 最 大 简约 (MP) 
树 ， 分 文 上 的 数字 表示 简约 性 分 析 1000 次 重复 得 到 的 靴 带 
分 析 支 持 率 (BS) ， 单 倍 型 的 地 理 分 布 情况 标 在 各 单 倍 型 后 


Fig.3 The most parsimonious (MP) tree based on eleven cpDNA 























haplotypes (H1-H11) of P. poissonii, numbers above the branches 
indicate MP bootstrap values (BS) based on 1000 replicates , 
and geographical distribution of each haplotype 


was marked out respectively 
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图 4 海 仙 报 春 两 个 叶绿体 DNA 序列 联合 分 析 的 11 个 单 倍 型 网 状 进 化 树 





圆圈 中 或 周围 的 数字 (Hl ~ H11) 表示 单 倍 型 





， 圆 点 表示 假定 的 单 倍 型 ， 每 一 分 支 表示 一 次 突变 








Fig.4 Nested cladogram of 11 haplotypes based on two cpDNA sequences of P. poissonii 





Cycles with numbers (H1-H11) denote haplotype. Dots re 


resent putative haplotypes. Each branch represents one mutation 


表 4 ”基于 对 海 仙 报 春 叶绿体 DNA 单 倍 型 数据 的 能 套 分 支 分 析 结果 ， 应 用 Templeton (2004) 
的 检索 表 分 析 得 到 的 一 系列 推论 
Table 4 Chain of inference from the nested clade analysis of the cpDNA haplotypes data of 


P. poissonii using Templeton’s (2004) inference key 








Clade Clade key ae Tierens 

1-4 1-19NO 异域 片段 化 Allopatric fragmentation 

1=5 Null hypothesis cannot be rejected Moving to nest clade 

2-4 1-19NO 异域 片段 化 Allopatric fragmentation 

2-5 1-19NO 异域 片段 化 Allopatric fragmentation 

3-1 1-19NO 异域 片段 化 Allopatric fragmentation 

3=3 1-19NO 异域 片段 化 Allopatric fragmentation 
Total clade 1-19NO 异域 片段 化 Allopatric fragmentation 





明 其 起 源 较 晚 ，Clade- 耳 的 4 个 单 倍 型 分 别 固定 
在 云南 中 部 的 昆明 和 景 东 居 和 群 及 四 川西 部 的 园 宁 
居 群 。 显 然 , 不 管 从 系统 树 的 结果 ,还 是 从 
NCA 分 析 结 果 看 来 ， 海 仙 报 春 形成 的 这 3 个 分 





























究 结果 表明 ， 在 海 仙 报 春分 布 范围 内 ， 居 群 间 的 
遗传 分 化 水 平 很 高 (Gs =0.916)， 同 样 ， 分 子 
变异 分 析 (AMOVA) 的 结果 表明 ， 居 和 群 间 的 遗 
传 变异 为 99. 08% , Fo, = 0. 99077 (P<0. 001), 











支 拥有 隔离 的 地 理 分 布 。 从 分 布 图 (图 1) 中 可 
以 看 出 ， 各 个 分 支 在 冰期 可 能 存在 各 自 的 避难 
所 ，NCA 结果 文 持 这 一 推断 : 隔离 分 化 是 主要 
的 分 此 模式 ， 即 在 各 分 支 存在 不 同 的 避难 所 。 
植物 居 群 空间 上 的 隔离 ， 突变、 环境 因子 造 
成 的 选择 差异 ,随机 遗传 漂 变 以 及 基因 流 的 隔离 
等 都 能 导致 居 群 遗传 结构 的 空间 异 质 性 ， 从 而 促 
进 居 群 的 分 化 (Jacquemyn 等 ,2004 ) 。 我 们 的 研 

















而 居 群 内 平均 遗传 多 样 性 H IE (0.078) 较 低 。 
可 以 看 出 海 仙 报 春 现 有 的 分 布 特征 为 : 居 群 间 具 
有 较 高 的 遗传 多 样 性 ， 而 居 群 内 的 遗传 多 样 性 较 
低 。 这 种 分 布 模式 与 大 多 数 温 带 植 物 在 间 冰 期 与 
冰期 的 扩张 现象 不 符 ， 后 者 是 分 享 同 一 个 避难 所 
或 来 自 同一 避难 所 的 居 和 群 表现 出 较 低 的 居 群 间 及 
居 群 内 遗传 多 样 性 ， 并 且 单 倍 型 组 成 较 单 一 
(Newton 等 , 1999; Petit 等 , 2003; Zhang 等 , 2005) 。 
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或 者 是 在 避难 所 的 居 群 和 冰期 后 扩散 的 居 群 的 遗 
传 多 样 性 不 同 ， 由 于 黄 基 者 效应 和 瓶颈 效应 使 冰 
期 后 扩散 形成 的 居 群 的 遗传 多 样 性 变 低 (Ibrahim 
等 , 1996; Avise, 2000; Hewitt, 1996, 2000) 。 显 
然 ， 在 我 们 的 研究 中 并 未 出 现 上 述 现象 ， 在 
NCA 分 析 推 新 形成 现 有 分 布 格局 的 主要 原因 是 
异域 片段 化 ， 即 各 居 群 间 的 基因 交流 受 限 。 海 仙 
报 春 的 基因 流 主 要 有 两 方面 ， 种 子 流 和 花粉 流 。 
种 子 缺乏 长 距离 传播 的 能 力 ， 而 花粉 流 主要 通过 
熊 蜂 等 昆虫 进行 传粉 。 叶 绿 体 DNA 在 海 仙 报 春 
中 为 种 子 传播 ， 居 群 间 的 平均 基因 流 为 0. 08, 
表明 不 同居 群 间 的 基因 交流 很 少 。 造 成 当前 这 种 
居 群 间 很 低 的 基因 流 和 很 高 的 遗传 分 化 现象 的 主 
要 原因 可 能 是 由 于 其 特殊 的 湿地 生境 ， 种 子 传播 
的 距离 相对 有 限 ， 而 花粉 传播 能 力 较 低 ， 加 之 居 
群 间 的 地 理 隔离 和 生境 片段 化 进一步 阻碍 了 居 群 
间 的 基因 交流 (Slatkin，1985; Koch 等 ，2006 ; 
Ehrich 等 , 2007) 。 

位 于 青藏 高 原 东南 部 的 横断 山地 区 的 植物 在 
第 四 纪 冰 期 时 ， 受 到 冰期 及 间 冰 期 的 影响 ,种群 
随 着 气候 的 波动 在 其 发 生地 反复 缩小 和 扩张 。 在 
冰期 时 由 于 该 区 域 没有 大 范围 的 冰 羡 ( 李 吉 均 
等 , 1996; 施 雅 风 , 2004)， 使 得 一 些 类 群 可 以 继 
续 在 其 发 生地 生存 ， 类 似 海 仙 报 春 的 高 山 类 群 可 
能 会 退缩 到 低 海 拔 地 区 渡 过 这 个 时 期 。 然 而 ， 在 
横断 山 多 数 地 区 都 存在 着 平行 走向 的 山脉 及 被 其 
分 隔 出 来 的 河谷 ， 对 于 生长 在 沼泽 或 湿地 的 海 仙 
报 春 ， 这 种 斑 块 状 的 生境 和 强烈 的 地 理 隔 离 限 制 
了 其 基因 交流 ， 加 之 海 仙 报 春 的 花粉 及 种 子 传播 
在 空间 距离 上 都 受 限 ， 花 粉 传播 距离 大 约 仅 有 
5.4 ~7.2m (Ishihama 等 , 2006) ， 显 然 难以 越过 
横 吾 的 高 山 。 所 以 ,不同 的 居 群 在 退缩 过 程 中 很 
难 融合 交流 形成 一 个 共同 的 避难 所 。 同 样 ， 居 和 群 
在 间 冰 期 也 很 难 在 扩张 的 过 程 中 进行 基因 交流 ， 
冰期 和 间 冰 期 的 不 断 反复 ， 导 致 种 群 反复 的 缩小 
和 扩张 ,但 是 它们 仍然 保留 着 各 自 相 对 隔离 的 地 
理 分 布 。 居 群 的 反复 缩小 和 扩张 使 单 倍 型 由 于 遗 
传 漂 变 而 随机 固定 在 不 同 的 居 群 中 ， 即 使 是 十 分 
临近 的 居 群 中 也 是 如 此 ， 遗 传 漂 变 降低 了 居 群 内 
的 遗传 多 样 性 ， 而 促进 了 居 群 间 的 遗传 多 样 性 。 
这 种 现象 导致 了 居 群 在 小 范围 内 扩张 和 集中， 从 
而 形成 了 现在 的 分 布 格局 ( Wang 等 ，2008; 
































































































































Meng 等 , 2008; Zhang 等 , 2005), 

对 海 仙 报 春 谱系 地 理学 的 研究 结果 表明 ， 该 
种 植物 在 冰期 可 能 存在 多 个 地 理 隔 离 的 避难 所 ， 
并 没有 表现 出 大 多 数 温带 植物 在 第 四 纪 冰 期 所 经 
历 过 的 种 群集 体 扩张 和 迁移 的 现象 。 由 于 横断 山 
地 区 特殊 的 、 复 杂 的 地 理 环 境 ， 如 平行 的 山脉 及 
山谷 形成 的 天 然 屏 障 ， 阻隔 了 种 群 的 基因 交流 及 
迁移 ， 形 成 了 海 仙 报 春 的 现代 分 布 。 本 研究 对 海 
仙 报 春 在 横断 山区 的 现代 分 布 格局 形成 的 历史 成 
因 进 行 了 讨论 ， 为 更 好 的 了 解 第 四 纪 冰 期 对 横断 
山地 区 植物 区 系 格局 的 形成 及 冰期 后 植物 居 群 的 
迁移 的 影响 ， 还 需 对 更 多 的 类 群 开展 研究 。 
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